function Q=findQ(S,J,I,M,p)
    Q=zeros(I,M)
        for m=1:M
        PIE=findPIE(J,m);
          for i=1:I
          %realPIE一个根据Sij删PIE某些行的程序
          PIE=realPIE(PIE,S,i);
          Q(i,m)=function_Q(i,m,J,S,p);
          end
        end
end